Quantum Computational Method of Finding the Ground State Energy and 

Expectation Values 
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We propose a new quantum computational way of obtaining a ground-state energy and expectation 
values of observables of interacting Hamiltonians. It is based on the combination of the adiabatic 
quantum evolution to project a ground state of a non-interacting Hamiltonian onto a ground state of 
an interacting Hamiltonian and the phase estimation algorithm to retrieve the ground-state energy. 
The expectation value of an observable for the ground state is obtained with the help of Hellmann- 
Feynman theorem. As an illustration of our method, we consider a displaced harmonic oscillator, 
a quartic anharmonic oscillator, and a potential scattering model. The results obtained by this 
method are in good agreement with the known results. 
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I. INTRODUCTION 

Quantum simulation might be a real application 
of medium-scale quantum computers with 50 — 100 
qubits As Feynman suggested, a quantum computer 
can simulate quantum systems better than a classical 
computer because it is also a quantum system Lloyd 
demonstrated that almost all quantum systems can be 
simulated on quantum computers Q . Abrams and Lloyd 
presented a quantum algorithm to find eigenvalues and 
eigenvectors of a unitary operator based on the quan- 
tum phase estimation algorithm Although it is an 
efficient quantum algorithm, there is room for improve- 
ment. First, one has to prepare an input state close to 
unknown eigenstates. Second, it has been little explored 
how to obtain physical properties except the energy spec- 
trum. 

In this paper, we propose a new refined quantum 
computational method to calculate the ground state en- 
ergy and expectation values of observables for interacting 
quantum systems. The main idea is as follows. Adia- 
batic turning on an interaction makes the ground state 
of a non-interacting system evolve to the ground state of 
an interacting system. During the adiabatic evolution, 
the phase estimation algorithm extracts the phase of an 
evolving quantum system continuously without the col- 
lapse of a quantum state. So the ground energy of an 
interacting system is obtained as a function of coupling 
strength. With the help of the Hellmann-Feynman the- 
orem Q, the expectation value of an observable for the 
ground state of an interacting system is obtained. As a 
test of our method, we simulate on classical computers 
three quantum systems: a displaced harmonic oscillator, 
a quartic anharmonic oscillator @ , and a potential scat- 



tering model 



II. METHOD 



Let us start with a brief review of Abrams and Lloyd's 
algorithm. Its goal is to find eigenvalues En and eigen- 
states \En) of a time-independent Schrodinger equation 

H\En) ^ En\En) . (1) 
Their key idea to solve ^ is to consider its time evolution 
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where \^ i) = X^n*^"!^") input or trial state. The 

information on eigenvalues En in the input state is trans- 
fered to index qubits by applying the quantum phase esti- 
mation algorithm. The measurement of the index qubits 
gives us a good approximation to En with probability 
|a„p, and makes j'l'/) collapse to \En)- It is instructive 
to compare ([2]) with the quantum Monte Carlo method 
which uses the imaginary time t = it to project the input 
state onto the ground state Q 
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First, in order to find the ground state energy, both ([2|) 
and dS]) require a good input state close to \Eo). If the 
input state does not contain the information about the 
ground state, both will fail. Second, for each run, while 
([2]) outputs En randomly, ([3]) produces Eq always. Fi- 
nally, is a real time evolution, however, ^ is the 
imaginary time evolution, i.e., a diffusion process, which 
is implemented by classical random walks. 

Our goal is to find a ground state energy with prob- 
ability 1 even if an input state contains little informa- 
tion on the ground state. Our method uses a real time 
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projection onto the ground state by adiabatically turn- 
ing on an interaction. Ortiz et al. suggested the use of 
the Gell-Mann-Low theorem to find the spectrum of a 
Hamiltonian with quantum computers 0, ml. Farhi et 
al. developed the adiabatic quantum computation [llj . 

We divide the Hamiltonian H into two parts: non- 
interacting Hamiltonian Hq and interaction Hi, H = 
Hq + Hi, As usual, it is assumed that the eigenvalues 
Wn and eigenstates \Wn) of Hq are known, Ho\Wn) = 
Wn\Wn)- We recast H to be time-dependent 



H{t) =Ho + f{t)Hi + E, 



(4) 



where a slowly varying function f{t) satisfies /(O) = 
and f{Tfi) = 1 with running time Tr. The role of the 
constant energy Ec will be explained later. As the inter- 
action is turned on slowly, the input state |Wo) evolves 
adiabatically to |i?o) 

fe-^^o Hindi' 1^^^ ^ Eo{t')dt' 1^^^ ^ 

where T is a time-ordering operator [l^. Notice the 
similarity and difference between and ([5]). The 

quantum phase estimation algorithm can extract the in- 
formation on Eq from ([5]). Since, during the adiabatic 
evolution, the quantum system is in an instantaneous 
ground state |£'o(i)) of H{t), one can apply frequently 
the phase estimation algorithm without the collapse of 
the quantum state to the exited states. 

Since the phase (p = Ent/h is defined in < < 2tt, 
the phase estimation algorithm gives us only the absolute 
value of an energy \En\. Its sign can be determined by 
adding Ec- When Eq is negative, while Wq is positive, 
Ec > Eo makes all the spectrum positive. Also Ec is 
useful for stabilizing the algorithm. If \Eq\ is close to 
zero, a long time needs to make the phase = Eot/h 
finite. 

The expectation value of an observable O can be ob- 
tained with the help of the Hellmann-Feynman theo- 
rem It states that if iJ(a)|£'„(a)) = E,i{a)\En{a)) 
with parameter a, then the following relation holds 



—— = iE^ia)\ — \E,,{a)) 



(6) 



By modifying the full Hamiltonian to have a linear cou- 
pling to O, H{t) = Ho + f{t){Hi + aO) + Ec, © becomes 



dEn{a) 



da 



= {E,,\0\E.„ 



(7) 



Therefore, the expectation value of an observable is 
obtained from a derivative of En{oL) at a = 0. In 
practice, ^ is obtained from a numerical approxi- 
mation {En\0\E„) = [En{a)- Er,{-a)]/2a + Oia'). 
This is comparable with an expectation estimation al- 
gorithm Notice that our scheme does not require 
the repeated measurements and the average over the in- 
dividual outcomes 



III. APPLICATION TO QUANTUM SYSTEMS 

A. Displaced harmonic oscillator 

As an illustration of our method, let us consider a sim- 
ple Hamiltonian, 



" 2m 



Hi = \x. 



(8) 



For convenience, we set h = in = lo = 1. It is well known 
that ([5]) is exactly solvable, a usual perturbation theory 
for it works well, and its ground state is a coherent state. 

The first step to quantum simulation is to map a phys- 
ical system to a qubit system. The position a; in Q is 
continuous, but qubits are discrete. A usual approach 
is to discretize x. Another way is to map the eigen- 
states |n) of Hq to the computational basis of qubits, 
|n) = |jjv-i) ® \jN-2) «) ■ • ■ |jo) with n = Jn-i 2^~^ + 

jN-2 2^"^ H hjo 2° and jk = or 1. Then Hq is given 

by a diagonal matrix. 
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The quantum dynamics of ^ was simulated on an NMR 
quantum computer by Somaroo et al. (l^ . On the other 
hand. Hi is written as a tridiagonal matrix, 

2"-2 

Hi^^Y. (v^i^)("+ii + v^i" + i)("i) • (10) 



n=0 



A quantum state \tp(t)) time t can be expressed in 

terms of \n), \i;{t)) = J2n=o "'n{t)\n). 

The adiabatic time evolution ^ is implemented by 
solving the time-dependent Schrodinger equation with 
the forth-order Runge-Kutta method on a classical com- 
puter. We take A^ = 3 ~ 6. We assume hat the phase es- 
timation algorithm is implemented very accurately. The 
adiabatic switching-on function f{t) used here is given by 
/(<) = i + itanh(20t/rfl.-10). One may expect it would 
take a long time for a quantum system to evolve adiabat- 
ically. However, in the case considered here, it takes the 
running time Tji = 15 Tq to obtain the ground state 
energy with accuracy AEq = lE^""^""^ - El™^\ < 10-^ 
where To = 2tt/llj is the period of the ground state of Hq, 
and E'n"™ is the numerical value. 
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Fig. [1] shows how the dynamical phase of the system 
changes as the interaction is slowly turned on. In Fig. [T] 

(a) , A = 0, and the oscillation period is To. Figs. [T] 

(b) and (c) show how the constant energy Ec is used 
to change the frequency corresponding to the ground 
state energy of an interacting Hamiltonian. In Fig. [T] 
(b), A = 0.9 and Ec = 0. So the frequency Eq/H become 
very low. However, in Fig. [1] (c), the constant energy 
Ec ~ 1/4 shifts the frequency so it can be easily mea- 
sured. In Figs.[l](d) and (e), we take A = V2 so the exact 
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FIG. 1: (color online). Re{ao(t)} as a function of t/To for 
(a) A = and £'c = 0, (b) A = 0.9 and = 0, (c) A = 0.9 
and = 0.25, (d) A = ^2 and E^ = 0, and (e) A = ^ 
and Ec = 1.0. Here TV = 4 and To = 2K:h/WQ = An with 
Wo = huj/2. 



energy is Eq = —1/2. Since the phase estimation algo- 
rithm produces only the absolute value of energy, \Eq\, 
the constant energy Ec is added in (H]) to decide its sign. 
In Fig. [1] (d), Ec = 0. At the end of running, the es- 
timated energy is Eq = 1/2. So the phase estimation 
algorithm fails to calculate the exact ground state en- 
ergy Eq = —1/2. However, in Fig. [T] (e), Ec = 1. the 
phase estimation algorithm gives us the energy 1/2. So 
we know that the exact energy is Eq = 1/2 — 1 = —1/2. 
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FIG. 2: (color online), (a) Probability p„(t) = |a„(t)P of 
qubits in the computational basis |n), and (b) the instanta- 
neous ground-state energy Eoit) in the unit of hu as a function 
of the dimensionless time I/Tr for A — y/Q. In (b), f{t) is an 
adiabatically switching-on function, (c) At t = Tr, pn obeys 
the Poisson distribution of a coherent state. 



FIG. 3: (color online). Ground state energy -Eo(q) in the unit 
of huj as a function of a for (a) A = and (b) A = 1. The 
points are numerical results. The red line in (a) is the plot of 
/(q) = i{a + l). In (b), the red one is the plot of g{a) = ^a. 



For any A, the ground state of ([5]) is a coherent state. 
As shown in Fig. [21 the probability that qubits are in the 
number state |7i) follows a Poisson distribution. So the 
ground state obtained by the quantum simulation might 
be called a pseudo-coherent state because it is defined on 
the truncated Hilbert space. It is a collective state of N 
qubits. 

The coherent state is also characterized by the min- 
imum uncertainties in x and p. Its mean square de- 



is 1/2 for any A. 



viation of x, Ax = (x ) — ( 
The ground state of ([8]) is displaced from the origin to 
x = —A. So (x) ~ —A. With the help of the Hellmann- 
Fcynman theorem, we calculate (x^) for A = and 
A = 1. To this end, the final Hamiltonian is modified 
as H{t) ^ Ho + f{t){Xx + ax^) -f Ec. Fig. [3] shows the 
ground state energy Eo{a) as a function of a. The deriva- 
tive of Eq (a) at a = gives us the expectation value of 

As illustrated in Fig. [31 we have 



X, [X 
„2 



dEpja) 
da 



a=0 



X') = 0.02/0.04 = 1/2 for A = 0. Thus Ax^ = 1/2. 
For A = 1, (^2) = 0.03/0.02 3/2. Again we have 
Aa;2 = 3/2 - 1 = 1/2. 
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B. Quartic anharmonic oscillator 

Let us consider an anharmonic oscillator, whose Hamil- 
tonian is given by 



Ho 



2m 



Hi = Ax* 



(11) 



where A > is the coupling constant. In their semi- 
nal paper Q , Bender and Wu showed that the Rayleigh- 
Schrodinger perturbation theory for (jlip becomes diver- 
gent for any A. Various non-perturbative methods have 
been applied to this simple model. 

One can write Hi = j{a^ + a)** = ^ + jV , where 



rich physics Q. The naive perturbation theory breaks 
down no matter small g is. For an attractive potential, 
i.e., g < 0, the lowest eigenstate of becomes a bound 
state. Also it exhibits the Anderson orthogonality catas- 
trophe [l^ which states that the ground state of Hq + Hi 
becomes orthogonal to the ground state of Hq in the ther- 
modynamic limit. 

We map the single-particle energy level of Hq to a 



computational basis. 



c}jC„|vac), where |vac) is a 
vacuum state. In p3)) . Hq can be written as a diago- 
nal matrix, (iJo)mn = ^n^mn- Whcrcas Hi are given by 
(i^l)m„ = <?/2^, which is more dense than (jlOp and ([12]). 



Vmn ^ Vin- ± l)(n ± 2)(n ± 3)(n + 2 ±2) 6 



+2(2n + 1 ± 2)^{n±l)in+l±l)6m,n±2 

+6n{n + 1) Sr,^,n ■ (12) 

The matrix of (fT^ is more denser than pH)) . So more 
qubits are used in (|12p in order to get the accurate energy. 




FIG. 4: (color online). For a quartic anharmonic oscillator, 
(a) Eo{X) in the unit of fioj and (b) Ax^ as a function of A. 
Here iV = 6 and = 15 To. 




Fig. m shows the ground state energy £'o('^) a-nd 
Ax^(A) as a function of A. For A = 2.0 and time step 
At = 5.0 X 10-^ we obtain Eq = 0.951568 472125, 
which is comparable to the best known results Eq ~ 
0.951568 472 722 For the calculation of Aa;^ wc 

obtain Eo{X, a) of Hq + XHi + ax'^ for a ^ ±0.001. Thus 
we have (x^) « [Eo{X, a) — Eq{X, —a) ] /2a. 

C. Potential scattering model 

Finally, we consider spinless electrons with a contact 
potential with Hamiltonian 

2" 

Ho = Y^ e,i cicn , Hi = cicn , (13) 

n— 1 n,m 

where e„ — (n — 1)A with level spacing A, cjj is a cre- 
ation operator, and g the coupling constant. Although 
this model is very simple and exactly solvable, it contains 



FIG. 5: (color online), (a) Energy levels En (in arbitrary 
unit) and fidelity Fn as a function of g. Here A'' = 6 and 
A = 10/64. 



As g is turned on adiabatically, the initial state |n) 
evolves to the final state \En{g)). We use the notation 
l-E'n(O)) = |e„) = \n). Fig. [5] (a) illustrates the single- 
particle levels En{g)- One see that there is one bound 
state with negative energy £'0(5) < for g < 0, but 
otherwise it is positive. Fig. [5] (b) shows the fidelity 
Fn{g) = \ {En{g)\'n)\'^ as a function of 9. Surprisingly, it is 
also calculated with the help of the Hellmann-Fcynman 
theorem. One can rewrite i^„((7) = (£'„(.g)|0|£'„(.g)) with 
O ~ I"-) ("-I- As shown in Fig. [5](b), the fidelity decrease 
more rapidly for g < than for g > 0. One can see that 
even single-particle levels for 5 = and g < become 
orthogonal. It is interesting that the fidelity between the 
interacting and non-interacting many-body ground states 
can be obtained from all the information of single-particle 
levels [13 ■ 
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IV. CONCLUSIONS 

In conclusion, we have proposed a new method to find 
the ground state energy by adiabatically turning on an 
interaction. The expectation values of an observable 
has been obtained by switching on a modified interac- 
tion which contains an observable and by applying the 
Hellmann-Feynman theorem. Our method has been suc- 
cessfully tested by solving three quantum systems. We 
expect that our method could be applied to the simula- 
tion of more interesting quantum systems. 

Finally, let us discuss the limits of our method. Our 
method is based on the combination of adiabatic quan- 
tum computation and the phase estimation algorithm. 
So, the computational resources needed to implement our 
method is approximately equal to the sum of those in- 
volved in adiabatic quantum computation and the phase 
estimation algorithm. The running time of the adiabatic 



evolution increases if the gap between the energy lev- 
els decreases. However, it is expected that the quantum 
Zeno effect [31 might release this limitation. A quantum 
state after applying a quantum phase estimation algo- 
rithm is approximately given by ~ ao|£'o)s |^^o)/ + 
ai\Ei)s where |aip = 1 — e and |aip = e for small e 
and subscripts "S" and "I" refer to the system and the in- 
dex qubits, respectively. The measurement on the index 
qubits gives us = |£'o)s l^^o)/ with high probabil- 
ity. The frequent applications of a quantum phase esti- 
mation algorithm and measurement on the index qubits 
could accelerate an adiabatic evolution. This will be in- 
vestigated in a future study. 
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